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I.  INTRODUCTION 


A  spinning  liquid  can  support  waves  within  its  interior.  These  waves 
have  been  studied  extensively  in  rotating  flows  and  are  commonly  known  as 
inertial  waves.1  For  a  cylinder  filled  with  a  viscous  liquid  that  is  spinning 
as  a  quasi-rigid  body,  the  wave  frequencies  (the  eigenfrequencies)  depend  upon 
the  cylinder  aspect  ratio  (half  hei ght/radius  =  c/a),  the  liquid  fill  ratio, 
and  Reynolds  number  (Re,  defined  in  subsequent  discussions).  When  a  liquid- 
filled  projectile  is  launched,  the  payload  must  be  spun  up  by  the  projec¬ 
tile.  The  circumferential  velocity  of  the  liquid  varies  during  spin-up,  and 
the  values  of  the  liquid  eigenfrequencies  will  also  depend  upon  this  time- 
dependent  velocity  distribution.  When  the  liquid  is  spinning  as  rigid  body, 
the  destabilizing  liquid  moment  can  be  computed  by  several  methods.2’9’4  All 
of  these  models  have  been  correlated  against  experimental  data  and  are  consis¬ 
tent  for  Re  >  1,000.  Detailed  information  about  spin-up  is  available  in  two 
areas:  unperturbed  velocity  histories5’6’7’8  and  eigenvalues  of  the  liquid 
during  spin-up.9  Recently,  however,  Murphy  has  developed  an  approximate 
scheme  for  the  determination  of  yaw  moments  during  spin-up.10 


1.  H,  P.  Greenspan,  The  Theory  of  Rotating  Fluids ,  Cambridge  University 
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Without  a  Critical  Layer,"  BRL  Technical  Report  in  preparation. 


This  report  describes  the  results  of  a  flight  test  where  spin-up  insta¬ 
bilities  were  obtained.  The  design  rationale  and  flight  data  are  discussed. 
The  design  of  the  test  was  heuristic  in  nature  since  at  that  time  a  predictive 
model  for  a  spin-up  liquid  moment  did  not  exist.  A  destabilizing  liquid 
moment  can  become  large  when  an  ei genf requency  (t  )  is  nearly  equal  to  the 

yaw  frequency  of  the  projectile  (t).  (The  mode  numbers  k  and  n  will  be 
explained  in  the  next  section.)  The  overall  stability  of  the  liquid-projec¬ 
tile  system,  however,  depends  upon  the  relative  size  of  the  destabilizing 
liquid  moment  and  the  aerodynamic  damping  moment. 


II.  SPIN-'JP  INSTABILITIES 
A.  Spin-up  Eigenfrequencies 

Assume  that  a  projectile  experiences  a  liquid-induced  moment  during  spin- 
up.  This  liquid  moment  may  be  short-lived  since  a  change  in  the  circumfer¬ 
ential  velocity  distribution  will  produce  a  change  in  the  liquid  eigenfre- 
quency.  If  the  total  spin-up  moment  is  sufficiently  large  to  induce  a  growth 
in  yaw,  then  the  growth  rate  would  have  the  following  functional  form: 

Yaw  Growth  Rate  =  F  (  c ,  Tkr)  ,  dtkn/dt)  (1) 


A  simple  analytic  expression  for  the  function  F  is  not  possible.  However,  it 
is  clear  that  the  yaw  growth  will  be  inversely  related  to  dikn/dt,  which  is 

the  time  variation  at  the  liquid  eigen  frequency  during  spin-up.  If  dtkn/dt  is 

small  the  spin-up  moment  will  be  large  since  the  resonant  behavior  will  per¬ 
sist  for  a  substantial  period  of  time.  On  the  other  hand,  if  dTkn/dt  is  large 

(this  is  normally  the  case  for  e3rly  flight  times)  then  the  time  duration  of 

the  moment  will  be  short  and  the  yaw  may  be  only  slightly  affected.  In  the 
design  of  actual  hardware,  the  physical  characteristics  of  the  projectile 
metal  parts  specify  x  and  can  not  be  varied  to  any  great  degree.  However,  the 
yaw  growth  rate  could  be  controlled  through  dr^/dt  which  is  a  sensitive 

function  of  Re  and  c/a.  Using  this  simple  concept,  spin-up  ei genf requency 
histories  were  computed  and  canister  geometries  and  liquids  were  chosen  to 
yield  particular  values  of  c/a  and  Re.  (Only  fill  ratios  of  100*  were  used.) 
Conditions  were  selected  which  would  produce  a  resonant  matching  between  a 
spin-up  ei genf requency  and  the  dimensionless  yaw  frequency  (fast  precession 
frequency/spin  frequency).  Many  combinations  of  Re  and  c/a  were  computed. 
Because  the  magnitude  of  the  liquid  moment  was  still  unknown,  a  prediction  of 
yaw  growth  could  not  be  made.  Using  the  method  outlined  in  Reference  9, 

histograms  of  spin-up  eigenfrequencies  were  computed.  A  typical  computation 

is  presented  in  Figure  1  for  a  155mm  projectile.  The  solid  line  is  a  computed 

spin-up  ei  genf  requency  history  for  c/a  =  4.972  and  Re  =  1,85  *  10^.  The 

computational  scheme  is  started  at  long  times  where  the  steady  state  eigen- 

frequency  is  known  from  any  of  the  linear  theories  and  is  advanced  backwards 

in  time.  At  some  very  early  time  (normally  less  than  1  second)  the  scheme 

will  not  converge,  but  this  minor  shortcoming  did  not  impact  the  design  of  the 
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experiment  (and  does  not  produce  any  practical  difficulties  with  the  general 
use  of  this  program  as  an  aid  to  projectile  designers). 

2 1 

Spin-up  is  controlled  by  the  Ekman  number  (E),  which  is  defined  as  v/c  $ 
and  represents  the  relative  magnitude  of  the  viscous  and  Coriolis  accelera¬ 
tions.  In  Reference  1  it  is  shown  that  the  characteristic  time  to  spin-up 
(the  e-folding  time)  is 

t(spin-up)  =  E”1'/2|r1  (2) 

Historically,  in  many  liquid-filled  shell  analyses,  a  Reynolds  number  (Re)  is 
2  *  2 

defined  as  a  hence,  Re  =  (c/a)  /E.  Strictly  speaking,  spiri-up  flows 

should  be  addressed  in  terms  of  E,  but  Re  will  be  used  for  historical  consis¬ 
tency.  The  liquids  used  in  the  flight  tests  were  silicone  oils  of  1  and  50 
centistokes  (cs)  viscosity.  Water  (under  standard  conditions)  has  a  kinematic 

2 

viscosity  of  approximately  1  cs  or  1  cm  /sec.  The  densities  of  the  1  and  50 
cs  oils  a 'e  0.853  and  0.960  gm/cc,  respectively. 

For  simplicity,  assume  that  the  spin  of  the  projectile  does  not  decrease 
from  the  launch  spin  rate.  Under  this  condition,  the  liquid  will  achieve  a 
state  of  rigid  body  rotation  at  sufficiently  long  times  and  the  liquid  eigen- 
frequencies  can  be  predicted  by  the  linear  theories.2  ’ 3  ,l4  The  response  of  the 
liquid  for  the  case  of  solid  body  rotation  is  often  resolved  into  modes.  (See 
Reference  3  for  a  detailed  discussion.)  This  modal  representation  simplifies 
the  description  of  the  three-dimensional  wave  motion  within  the  liquid.  The 
radial  mode  number  is  n(l,2,3...),  and  the  longitudinal  mode  number  is 
k(l,3,5...).  In  order  to  produce  a  destabilizing  yaw  moment,  the  azimuthal 
mode  number  must  be  unity;  the  azimuthal  mode  number  is  implicitly  understood 
to  always  be  unity.  As  with  most  harmonic  oscillations,  only  the  fundamental 
modes  (low  mode  numbers)  are  of  practical  importance. (Higher  modes  are  usually 
damped  by  viscous  effects.)  Unfortunately,  such  a  judgment  is  based  upon  ex¬ 
perience  or  intuition  for  the  particular  problem  at  hand.  For  example,  the 
computed  data  in  Figure  1  were  based  upon  k =5  and  n=l.  It  is  convenient  to 
identify  a  spin-up  ei genf requency  history  by  the  mode  numbers  (k,n)  of  the 
steady  state  response. 

B.  Selection  of  Hardware 


A  single  occurrence  of  a  spin-up  instability  was  documented  in  Reference 
11.  Only  a  single  configuration  (a  modified  M678  with  a  half-caliber  base, 
c/a  =  4.972,  fill  ratio  =  100%,  and  Re  =  1.86  x  106  )  was  tested.  Results 
from  a  three  round  group  indicated  only  one  flight  instability;  however,  the 
resonance  time  as  predicted  by  the  spin-up  ei genf requency  history  correlated 
well  with  the  yawsonde  data.  The  spin-up  ei genf requency  history  of  this 


12.  P.P.  D1  Amico,  W.  Clay,  and  A.  Mark,  "Yauisonde  Data  for  MA37-Type  Projec¬ 
tiles  With  Application  to  Rapid  Spin  Decay  and  Stewartson-i'ype  Spin-up 
Instabilities,"  BRL  Memorandum  Report  ARB  RL  -MR-  03027,  ,'ur.e  1980 
(AD  A089646). 
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configuration  (called  an  M687L)  is  shown  in  Figure  1  and  was  the  basis  of  the 
selection  of  hardware  for  this  present  test.  The  (1,5)  mode  of  Figure  1  shows 
a  monotonical ly  decreasing  behavior  with  time.  (Actually,  for  time  equal  to 
zero,  the  ei genf requency  is  zero  and  a  rapid  increase  in  ei genf requency  occurs 
for  very  early  times.)  If  the  ratio  of  the  fast  precessional  frequency  to  the 

spin  frequency  (<t^/|)  is  equal  to  the  spin-up  ei  genf  requency ,  then  resonance 
occurs  and  a  growth  in  yaw  is  possible  (but  is  not  guaranteed).  Typically,  an 
M687  projectile  launched  at  transonic  velocities  will  have  ^/|  =  0.08. 

Hence,  from  Figure  1  resonance  should  occur  for  t  =  at  about  7  seconds. 
In  order  to  produce  a  larger  moment  by  forcing  dr^  /dt  to  be  smaller,  the 

time  of  the  resonant  matching  should  be  larger  than  7  seconds.  This  could  be 
accomplished  by  increasing  Re  or  c/a.  An  increase  in  either  of  these  vari¬ 
ables  would  increase  the  spin-up  time  and  result  in  both  a  longer  time  to 
resonance  and  a  small  r  value  of  dtkn  /dt.  From  a  practical  standpoint,  it 

would  be  advantageous  to  change  c/a  and  then  test  several  Re  values.  This 
approach  was  taken  for  the  (1,5)  mode.  It  was  also  possible  to  investigate 
the  (1,3)  mode  by  reducing  the  canister  length.  The  (1,3)  mode  will  produce 
the  same  liquid  frequency  as  the  (1,5)  mode,  but  the  effect  of  reduced  liquid 
mass  and  lower  longitudinal  mode  number  will  be  observed.  Figures  2  and  3 
show  ei genf requency  histories  for  (1,5)  and  (1,3)  modes  for  two  Reynolds 
numbers  (1  and  50  cs  oil). 

Hardware  was  constructed  for  the  aspect  ratios  shown  in  Figures  2  and 
3.  A  sectioned  view  of  a  standard  M687  is  shown  in  Figure  4.  Single  piece 
canisters  were  fabricated  utilizing  M687-type  components.  The  (1,3)  canisters 
were  constructed  with  the  center  of  the  cavity  in  the  middle  of  the  payload 

section  with  very  thick  endwalls.  Bases  with  quarter-caliber  boattails 

(typical  of  the  standard  M483A1  family  of  shell)  were  utilized,  rather  than 
the  half-caliber  boattails  that  were  used  for  the  M687L  shell.  Table  1  lists 
the  projectile  types  and  pertinent  physical  characteristics. 


TABLE  1.  PROJECTILE /CANISTER  DESCRIPTIONS 


Projectile 

Type 

Year 

Boattai 1 

(cal ) 

c/a 

a 

(cm) 

m^a2* 
(kg»m2  ) 

Fill  Ratio 

(%) 

M687 

1974 

1/4 

4.390 

5.37 

0.012 

80 

XM687L 

1978 

1/2 

4.973 

5.37 

0.014 

100 

XM678A 

1982 

1/4 

5.200 

5.14 

0.012 

100 

XM687B 

1982 

1/4 

3.120 

5.14 

0.007 

100 

*Dased  upon 

the  density 

of  water 

and  for  a 

fill  ratio 

of  100%. 
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III.  RESULTS 


A.  Test  Site 


A  test  program  was  conducted  at  NASA  Wallops  Island  (WI),  Virginia,  during 
September  1982.  A  ground  receiving  station  was  operated  by  BRL  on  the  Pad  2 
launch  site.  Radar  tracking  and  backup  telemetry  coverage  were  provided  by 
Wl  personnel.  An  M185  tube  with  a  standard  muzzle  break  was  used  in  a  univer¬ 
sal  mount.  Launch  velocities  were  obtained  by  smear  cameras.  Table  2  gives  a 
round-by-round  summary  of  the  test  program.  All  projectiles  were  fitted  with 
fuze-configured  yawsondes. 12  Yawsonde  data  are  in  the  form  of  solar  aspect 
angle  (Sigma  N,  which  is  related  to  the  peak-to-peak  yawing  motion  about  the 
trajectory)  and  time  rate  of  change  of  the  Eulerian  roll  angle,  $  (when  the 
yaw  is  small  $  is  nearly  the  spin;  the  data  are  simply  labeled  spin).  The 
yawing  motion  of  the  projectile  consists  of  two  modes  of  precession,  fast  and 
slow.  Typical  values  for  the  ^a^t  and  slow  modes  (at  transonic  launch  veloci¬ 
ties  for  155mm  shell)  are  10  and  1  Hz,  respectively.  The  amplitude  of  the 
fast  precessional  mode  is  destabilized  by  the  liquid  payload. 

TA8LE  2.  ROUND-BY-ROUND  SUMMARY3 


Date 

Round 

Type 

BRL 

Number 

WI 

Number 

Muzzle 

Velocity 

(m/s) 

Projecti  Te 
Mass 
(k<j) 

Range 

(m) 

FMAb 

(deg) 

Resonant0 

Time 

(sec) 

23 

Sep 

82 

Al-1 

1649 

El-397 

352.7 

35.6 

8,549 

— 

23 

Sep 

82 

A 1  -2 

1765 

El -398 

351.7 

35.7 

8,686 

4 

3.0 

23 

Sep 

82 

Al-3 

1775 

El-399 

353.6 

35.7 

8,961 

3.5 

--- 

29 

Sep 

82 

A 1-4 

1808 

El- 402 

352.3 

35.8 

9,052 

1.5 

5.0 

29 

Sep 

82 

A50-1 

1810 

El-403 

348.7 

36.4 

8,793 

1.5 

5.0 

29 

Sep 

82 

B 1  —  1 

1814 

El-405 

353.3 

34.7 

8,779 

1.0 

4.0 

29 

Sep 

82 

Bl-2 

1815 

El-406 

354.5 

34.8 

8,671 

1.0 

3.5 

29 

Sep 

82 

Bl-3 

1816 

El-407 

353.9 

34.7 

8,595 

1.5 

4.0 

29 

Sep 

82 

B50-1 

1848 

El-408 

352.3 

35.3 

8,869 

— 

3  All  rounus  were  launched  at  a  quadrant  elevation  of  750  mils  with  charge  4W. 
b  FMA  is  defined  as  half  of  the  first  maximum  amplitude  observed  within  the 
yawsonde  data. 

c  Time  at  which  growth  of  the  fast  mode  precession  was  observed  within  the 
yawsonde  data. 


12.  W.  H.  Mer-nagen  ana  U.  H.  Clay,  "The  Design  of  a  Second  Generation  Yaw- 
sonde,"  SRL  Memorandum  Report  ARBRL-MR-2368 ,  April  1974  (AD  7800641. 
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R. 


Yawsonde  Data 


Table  3  lists  the  number  of  unstable  rounds  within  each  group.  (Also  listed 
are  Reynolds  numbers  and  viscosities.)  The  data  will  be  discussed  by  projec¬ 
tile  type:  A  and  B. 


TABLE  3.  NASA  WALLOPS  ISLAND  TEST  --  SEPTEMBER  1982 


Projectile  Type 

Vi scosity 
(cs) 

Re 

Rounds  Tested 
Unstable/Total 

XM687A 

1 

1.99  x  106 

3/4 

XM687A 

50 

3.98  x  104 

1/1 

XM687B 

1 

1.99  x  106 

3/3 

XM678B 

50 

3.98  x  104 

0/1 

1.  Type  A  Projectiles 

Figures  5a  and  5b  show  limited  yawsonde  data  for  Round  Al-1.  The 
telemetry  data  contain  spurious  pulses  of  noise  producing  erroneous  yaw  data. 
(Spin  data  were  not  sensitive  to  the  noise.)  It  is  clear,  however,  that  a 
small  launch  yaw  existed  and  larger  yaw  resulted  downrange.  No  large  reduc¬ 
tion  in  spin  was  observed.  Figures  6a,  6b,  and  6c  give  the  yawsonde  data  for 
Round  Al-2.  Upon  launch  the  yaw  was  dominated  by  the  slow  precessional 
mode.  The  amplitude  of  the  fast  precessional  mode  decayed  until  approximately 
3  seconds  where  growth  was  seen.  The  slow  mode  then  damped,  while  the  fast 
mode  slowly  grew  to  an  amplitude  of  5  degrees  at  15  seconds.  During  the 

remainder  of  the  flight,  the  fast  mode  continued  to  grow,  but  catastrophic 

levels  of  yaw  were  not  observed.  The  spin  data  show  a  modulation  (charac¬ 
teristic  of  yawing  projectiles),  but  large  despin  did  not  occur.  Figures  7a, 
7b,  and  7c  show  data  for  Round  Al-3.  This  round  was  stable.  However,  within 
Figure  7b  the  fast  precessional  mode  damps  between  0-2.5  seconds  and  then 
grows  slightly  in  the  2.5-15  second  time  frame.  The  fast  mode  damped  during 
the  down  leg  of  the  trajectory  unlike  the  previous  two  rounds.  No  other 

anomalies  were  noted.  The  data  for  Round  Al-4  are  given  in  Figures  8a,  8b, 
and  8c.  These  data  are  essentially  identical  to  Round  Al-2.  Very  early  data 
indicate  a  small  FMA  of  1.5  degrees  that  is  dominated  by  the  fast  mode 
(Typically,  the  early  yaw  contains  both  slow  and  fast  modes  and  is  slightly 

dominated  by  the  slow  mode.)  The  fast  mode  damped  subsequent  launch  and  began 
to  grow  at  approximately  5  seconds.  Figures  9a,  9b,  and  9c  show  data  for  the 
only  50  cs  projectile,  Round  A50-1.  The  data  were  very  similar  to  the 
previous  unstable  rounds  which  were  loaded  with  1  cs  oil.  The  very  early  yaw 
was  dominated  by  the  fast  mode,  which  rapidly  damped  and  then  grew  beginning 
at  about  5.0  seconds.  No  unusual  spin  data  were  observed. 
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2.  Type  B  Projectiles 


Figures  10a,  10b,  and  10c  show  the  yawsonde  data  for  Round  Bl-1. 
Early  data  show  an  FMA  of  1.0  degree  with  rapid  damping  for  the  fast  pre- 
cessional  mode.  The  fast  mode  then  dramatically  grew  at  approximately  4.0 
seconds  and  resulted  in  a  yaw  amplitude  of  5  degrees.  The  spin  data  were 
similiar  to  the  A  series  shell.  Figures  11a,  lib,  and  11c  give  yawsonde  data 
for  Round  Bl-2.  The  initial  yawing  motion  was  small  (FMA  =  1.0  degree  and 

dominated  by  the  slow  mode),  but  the  fast  precession  began  to  increase  at  3.5 

seconds.  These  data  were  very  similar  to  Round  Bl-1.  Data  for  Round  Bl-3  are 
given  in  Figures  12a,  12b,  and  12c.  An  FMA  of  1.5  degrees  was  observed,  and 
the  fast  precessional  mode  began  to  increase  in  amplitude  at  approximately  4.0 
seconds.  The  spin  data  were  typical  of  the  previous  rounds.  Data  for  the  50 

cs  projectile  in  this  series  are  in  Figures  13a  and  13b.  No  unusual  yaw  or 

spin  behavior  was  noted. 


IV.  DISCUSSION 


A  short  description  of  the  spin-up  analysis  by  Murphy10  will  be  presented 
prior  to  a  discussion  of  the  data.  The  equations  governing  the  liquid  eigen- 
frequencies  during  spin-up  indicate  the  presence  of  a  critical  layer  within 
the  flow  domain.  This  critical  layer  is  a  region  where  the  local  oscillation 
frequency  (during  the  spin-up  process)  matches  the  rotational  velocity  of  the 
main  azimuthal  flow.  Considerable  effort  must  be  expended  to  properly  treat 
the  critical  layer,  but  Murphy  has  suggested  a  model  that  is  applicable  when 
the  critical  layer  is  not  a  dominant  factor.  This  model  can  be  used  to 

compute  yaw  moments  during  spin-up  except  at  early  times.  (This  is  usually 
where  the  critical  layer  is  important.)  Flight  data  for  the  50  cs  rounds  can 
be  treated  by  this  theory,  while  the  higher  Re  cases  will  require  a  more 
complete  theory.  It  must  also  be  stated  that  the  present  spin-up  analysis  has 
been  linearised  for  small  amplitudes  and  is  a  quasi-steady  state  approach  that 
does  not  reflect  unsteady  effects  such  as  the  coning  history,  etc.  Qualita¬ 
tive  comparisons  between  this  spin-up  theory  and  the  flight  data  should  be 

made  as  a  guide  to  the  development  of  more  complete  analyses  and  as  a  basis 

for  future  flight  tests. 


Assuming  only  the  fast  yaw  mode  is  important,  the  liquid  moment  is 
assumed  to  have  the  following  form. 


Liquid  Moment 


m,  a4$ 


23,2 


T1  K1 


U 
e  1 


(C 


LSM, 


+  i  C 


LIM, 


(3) 


Both  moment  coefficients  can  be  computed  during  spin-up  by  the  methods  within 
Reference  10,  but  the  projectile  yaw  Is  controlled  only  by  CL$M  .  Computed 

results  are  shown  in  Figures  14a  and  14b  for  the  A-50  and  B-50  cases  (for  c/a 
-  5.20  and  c/a  -  3.12  at  Re  =  3.98  x  104 ) .  The  independent  variable  shown  is 
$  =  4>t ,  the  roll  angle  from  time  zero.  is  a  function  of  the  yaw  growth 


'1 


rate  and  the  coning  frequency.  These  calculations  for  are  for  t  = 

0.090.  The  solid  lines  are  the  total  moment  coefficient  due  to  both  pressure 
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and  viscous  effects,  while  the  dashed  lines  indicate  only  the  part  due  to 
pressure.  The  area  under  the  curve  in  Figures  14a  or  14b  represents  the 
angular  impulse  per  radian  or 


Liquid 

IT 


Moment  \ 


1 


dt  =  mj  a25>2 


d$ 


(4) 


In  a  qualitative  sense,  the  flight  data  suggest  that  the  unstable  A  series 
should  have  a  larger  angular  impulse  than  the  stable  B  series.  The  angular 
impulse  should  be  computed  over  the  same  time  domain  for  such  a  comparison; 
however,  this  was  not  possible.  The  presence  of  the  critical  layer  prevented 
proper  estimations  for  C|_sM^  over  the  same  times  for  both  aspect  ratios. 


Hence,  a  simpler  criterion  was  sought, 
case  only  the  angular  impulse  due  to 


It  is  assumed  that  in 
the  local  maxima  in 


the 

C 


actual 


LSM 


1 


Is 


flight 

impor¬ 


tant.  Tnis  could  be  justified  somewhat  since  aerodynamic  damping  has  been 
omitted,  but  admittedly  a  complete  angular  impulse  curve  would  be  prefer¬ 
able.  A  simple  criterion  can  be  established  by  the  use  of  a  so-called 
"quality  factor"  or  "Q,"  which  is  quite  common  in  the  analysis  of  circuit 
resonances  in  electrical  engineering.  The  bandwidth  of  the  resonance  is 
determined  from  0.707  times  the  maximum  response.  In  the  problem  at  hand,  the 
Phi  intercepts  for  0.707  (CLSM  )  max  are  easily  t0und.  The  bandwidths,  B,  are 

shown  on  Figures  14a  and  14b.  The  areas  under  each  curve  (as  defined  by  the 
bandwidth)  were  measured  and  were  found  to  be  approximately  equal.  The  actual 


moment  then  scales  as  m^a2$2  .  For  c/a  =  5.20  the  total  moment  experienced 

during  spin-up  would  be  largor  than  for  c/a  =  3.12,  since  all  quantities  are 
the  same  except  for  mL. 


A  summary  of  computed  and  observed  data  is  presented  in  Table  4.  The 
computed  resonance  time  was  selected  at  the  spin-up  eigenfrequency  of  0.091. 

This  value  of  4>j/$  (fast  precessional  frequency/spin  frequency  -  10  Hz/110  Hz) 

was  taken  as  a  nominal  value  for  all  rounds. 


Qualitative  agreement  was  found  between  the  computed  times  for  resonance 
and  the  observed  times  for  yaw  growth  for  the  first  three  cases.  For  the  B-50 
cs  case  where  dTkn  /dt  was  the  steepest,  yaw  growth  was  not  observed,  which  Is 

qualitatively  correct.  It  is  interesting  to  note  that  the  reduction  in  the 
liquid  mass  for  the  B  series  did  not  yield  stable  flights  for  the  higher 
Reynolds  number  case.  However,  the  reduction  in  the  liquid  mass  and  a  larger 
value  of  dtkn  /dt  (caused  by  a  lower  Re  value)  produced  a  stable  flight  for 

the  single  B-50  cs  round.  This  is  In  qualitative  agreement  with  Murphy's 
spin-up  moment  calculation.  In  general,  the  observed  unstable  yaw  behavior 
was  not  catastrophic  nor  were  range  losses  extreme.  However,  the  yaw  that  was 
experienced  is  unacceptable  In  that  the  drag  is  increased  compared  to  the  drag 
of  solid  payload  rounds  with  the  same  exterior  shape.  It  is  possible  that 
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a  slight  reduction  in  the  aerodynamic  damping  (perhaps  stimulated  by  a  lower 
atmospheric  density  or  a  half  caliber  boattail)  would  have  produced  large  yaw 
and  substantial  range  losses.  Although  the  projectile  sample  sizes  were  very 
small,  it  was  encouraging  that  the  projectile  flight  behavior  (stable  or 
unstable)  was  uniform. 


TABLE  4.  COMPARISON  OF  COMPUTED  AND  OBSERVED  DATA 


Aspect 

Ratio 

Round 

Type 

Computed 

Time 

(Sec) 

Resonant 

Time 

(Sec) 

Reynolds 

Number 

dt.  /dt 
kn 

Number  of  Rounds 
Unstable/Total 

5.200 

A-l 

2-3 

3-5 

1.99  *  106 

-0.008 

3/4 

5.200 

A- 50 

2-3 

5 

3.98  x  106 

-0.44 

1/1 

3.200 

6-1 

5-6 

3.5-4 

1.99  x  106 

-0.013 

3/3 

3.200 

B-50 

1-2 

3.98  x  106 

-0.66 

0/1 

V.  CONCLUSIONS 

A  liquid- filled  projectile  can  become  unstable  due  to  a  resonance  between 
the  fast  precessional  frequency  and  a  liquid  eigenfrequency.  During  the 
liquid  spin-up  process  when  the  eigenfrequency  varies  slowly,  the  destabiliz¬ 
ing  moment  can  become  large.  Conversely,  If  the  eigenfrequency  varies  rapidly 
with  time,  then  resonance  effects  are  reduced  and  destabilizing  moments  may  be 
avoided.  With  the  use  of  a  spin-up  eigenfrequency  code,  projectile  hardware 
was  designed  and  tested  to  Investigate  these  effects.  Flight  data  were  In 
qualitative  agreement  with  the  eigenfrequency  estimates  and  Indicated  that  the 
magnitude  of  spin-up  moments  can  be  reduced  by  a  rapid  time  variation  of  the 
liquid  eigenfrequency  in  the  vicinity  of  resonance.  An  approximate  spin-up 
theory  was  applied  to  two  flight  cases  and  was  also  In  qualitative  agreement. 
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Figure  8c.  Spin  Versus  Time  for  Al-4 
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Figure  12a.  Sigma  N  Ve  ^  Time  for  Bl-3  (0-40  Sec). 
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Figure  14a.  C.  <-M  During  Spin-Up  for  c/a 
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Figure  14b.  C,  During  Spin-Up  for  c/a  ■  3.12. 
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Armament  R&D  Center 
US  Army  AMCCOM 
ATTN:  DRSMC-LCB-TL(D) 
Watervliet,  NY  12189 

1  Commander 

US  Army  Aviation  Research 
and  Development  Command 
ATTN:  DRDAV-E 
4300  Goodfellow  Blvd 
St.  Louis,  MO  63120 

1  Director 

US  Army  Air  Mobility  Research 
and  Development  Laboratory 
Ames  Research  Center 
Moffett  Field,  CA  94035 

1  Commander 

US  Army  Communications  Research 
and  Development  Command 
ATTN:  DRSEL-ATDD 
Fort  Monmouth,  NJ  07703 

1  Commander 

US  Army  Electronics  Research 
and  Development  Command 
Technical  Support  Activity 
ATTN:  DELSD-L 
Fort  Monmouth,  NJ  07703 

2  Commandant 

US  Army  Infantry  School 
ATTN:  ATSH-CD-CSO-OR 
Fort  Benning,  6A  31905 

1  AFWL/SUL 

Ki rt land  AFB,  NM  87117 
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DISTRIBUTION  LIST 


No.  Of 

Copies  Organi zation 

1  Commander 

US  Army  Missile  Command 

ATTN:  DRSMI-R 

Redstone  Arsenal ,  AL  35898 

1  Commander 

US  Army  Missile  Command 
ATTN:  DRSMI-YDL 
Redstone  Arsenal ,  AL  35898 

1  Commander 

US  Army  Missile  Command 
ATTN:  DRSMI-RDK,  Mr.  R.  Deep 
Redstone  Arsenal ,  AL  35898 

1  Commander 

US  Army  Tank  Automotive 
Command 

ATTN:  DRSTA-TSL 
Warren,  MI  48090 

1  Pi  rector 

US  Army  TRADOC  Systems 
Analysis  Activity 
ATTN:  ATAA-SL 
White  Sands  Missile  Range 
NM  88002 

1  Commander 

US  Army  Dugway  Proving  Ground 
ATTN:  STEDP-MT 

Mr.  G.  C.  Travers 
Dugway,  UT  84022 

I  Commander 

US  Army  Yuma  Proving  Ground 
ATTN:  STEYP-MTW 
Yuma,  AZ  85365 

3  Direct- or 

Sandia  National  Laboratory 
ATTN:  Mr.  H.  R.  Vaughn 
Dr.  W.  Oberkampf 
Mr.  F.  G.  Blottner 
Albuquerque,  NM  87115 


No.  of 

Copies  Organization 

1  Commander 

Naval  Surface  Weapons  Center 
ATTN:  Dr.  W.  Yanta 

Aerodynamics  Branch 
K-24,  Building  402-12 
White  Oak  Laboratory 
Silver  Spring,  MD  20910 

1  AFATL  (DLDL,  Dr.  D.C.  Daniel) 
Eglin  AFB,  FL  32542 

1  Aerospace  Corporation 
Aero-Engineering  Subdivision 
ATTN:  Walter  F.  Reddall 

El  Segundo,  CA  90245 

2  Director 

National  Aeronautics  and  Space 
Administration 
Ames  Research  Center 
ATTN:  Dr.  P.  Kutler 
Dr.  T.  Steger 
Moffett  Field,  CA  94035 

1  Director 

National  Aeronautics  and  Space 
Administration 
Langley  Research  Center 
ATTN:  Tech  Library 
Langley  Station 
Hampton,  VA  23365 

1  Director 

National  Aeronautics  and  Space 
Administration 

Marshall  Space  Flight  Center 
ATTN:  Dr.  W.  W.  Fowl  is 
Huntsville,  AL  35812 

2  Cal  span  Corporation 
ATTN:  G.  Homicz 

W.  Rae 
P.0.  Box  400 
Buffalo,  NY  14225 
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DISTRIBUTION  LIST 


No.  of 

Copies  Organization 

1  Martin-Marietta  Corporation 
ATTN:  S.H.  Maslen 
1450  S.  Rolling  Road 
Baltimore,  MD  21227 

?  Rockwell  International 
Science  Center 
ATTN:  Dr.  V.  Shankar 

Dr.  S.  Chakravarthy 
1049  Camino  Dos  Rios 
Thousand  Oaks,  CA  9136C 

1  Arizona  State  University 
Department  of  Mechanical  and 
Energy  Systems  Engineering 
ATTN:  G.P.  Neitzel 
Tempe,  AZ  85281 

1  Massachusetts  Institute  of 
Technology 
ATTN:  H.  Greenspan 
77  Massachusetts  Avenue 
Cambridge,  MA  02139 

1  North  Carolina  State  University 
Mechanical  and  Aerospace 
Engineering  Department 
ATTN:  F.F.  DeJarnette 
Raleigh,  NC  27650 

1  Northwestern  University 
Department  of  Engineering 
Science  and  Applied 
Mathematics 
ATTN:  Dr.  S.H.  Davis 
Evanston,  II  60201 

1  Notre  Dame  University 
Department  of  Aero  Engr 
ATTN:  T.J.  Mueller 
South  Bend,  IN  46556 

1  Rensselaer  Polytechnic 
Institute 

Department  of  Math  Sciences 
ATTN:  R.C.  Diprima 
Troy,  NY  12181 


No.  of 

Copies  Organization 

.  .X  University  of  California  -  • 

Davis 

ATTN:  H.A.  Dwye- 
Davis,  CA  95616 

1  University  of  Colorado 
Department  of  Astro-Geophysics 
ATTN:  E.R.  Benton 
Boulder,  CO  80304 

2  University  of  Maryland 
ATTN:  W.  Melnik 

J.D.  Anderson 
College  Park,  MD  20742 

1  University  of  Maryland  - 
Baltimore  County 
Department  of  Mathematics 
ATTN:  Dr.  Y.M.  Lynn 
5401  Wllkens  Avenue 
Baltimore,  MD  21228 

1  University  of  Santa  Clara 
Department  of  Physics 
ATTN:  R.  Greeley 
Santa  Clara,  CA  95053 

2  University  of  Southern 

California 

Department  of  Aerospace 
Engineering 
ATTN:  T.  Maxworthy 
P.  Weidman 

Los  Angeles,  CA  90007 

2  University  of  Rochester 
Department  of  Mechanical 
and  Aerospace  Sciences 
ATTN:  R.  Gans 

A.  Clark,  Jr. 

Rochester,  NY  14627 

1  University  of  Tennessee 
Department  of  Physics 
ATTN:  Prof.  W.E.  Scott 
Knoxville,  TN  37916 
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DISTRIBUTION  LIST 


No.  of 
Copies 


Organization 


3_  Virginia  Polytechnic  Institute 
"  and  State  University  .  .1 

Department  of  Aerospace 
.  Engineering 

ATTN:  .  Tech  Library  _____  __ 
"  "  Dr.  W.  Sarlc  : 

Dr.  T.  Herbert 
Blacksburg,  VA  24061 


1  University  of  Wyoming 
ATTN:  D.L.  Boyer 
University  Station 
Laramie,  WY  82071 


Aberdeen  Proving  Ground 

Director,  USAMSAA 
ATTN:  DRXSY-D 

DRXSY-MP,  Mr.  H.  Cohen 

Commander,  USATECOM 
'  ATTN:  DRSTE-TO-F 
•  DRSTE-CM-F 

Mr.  Gibson  (2  cys) 


PM-SMOKE  . 

ATTN:  DRCPM-SMK-M 

Cdr,  CRDC,  AMCCOM 
ATTN:  DRSMC-CLN 

Mr.  w.  Dee 
Mr.  C.  Hughes 

.  Mr.  F.  Dagostin 

Mr.  C.  Jeffers 
Mr.  L.  Shaft 
Mr,  M.  Parker 
DRSMC-CL9-PA 

Mr.  M.  C.  Miller 
DRSMC-CLJ-IL 
DRSMC-CLJ-L 
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Please  take  a  few  minutes  to  answer  the  questions  below;  tear  out 
this  sheet,  fold  as  indicated,  staple  or  tape  closed,  and  place 
in  the  mail.  Your  comments  will  provide  us  with  information  for 
improving  future  reports. 

1 .  BRL  Report  Number _ '  _ 

2.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related 
project,  or  other  area  of  interest  for  which  report  will  be  used.) 


3.  How,  specifically,  is  the  report  being  used?  (Information 
source,  design  data  or  procedure,  management  procedure,  source  of 
ideas,  etc.)  _  _ _ 


4.  Has  the  information  in  this  report  led  to  any  quantitative 
savings  as  far  as  man-hours /contract  dollars  saved,  operating  costs 
avoided,  efficiencies  achieved,  etc.?  If  so,  please  elaborate. 


nts  (I 
and  fu 
re  usa 


<>.  If  you  would  like  to  be  contacted  by  the  personnel  who  prepared 
this  report  to  raise  specific  questions  or  discuss  the  topic, 
please  fill  in  the  following  information. 

Name: 


Telephone  Number: 
Organization  Address: 


